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Abstract 

The weight method of implementing the BE effect into Monte Carlo 
generators is discussed and presented in some detail. We show how the 
choice of free parameters and the definition of "direct" pions influence 
the results for the hadronic decays. 
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1 Introductory remarks 



The effect of Bose-Einstein symmetrization (BE effect) in the two-particle 
correlation spectra depends on the shape and size of the source. This allowed 
to estimate the source parameters of astronomical sources via the so-called 
Hanbury-Brown and Twiss effect The analogous estimates in particle 
physics are much more involved H, |]. In fact, the applicability of the 
standard analysis assuming incoherent production in particle collisions was 
questioned recently and an alternative approach was presented ||^. The im- 
plementation of BE effect into Monte Carlo generators modelling multiple 
production is particularly difficult, as the symmetrization should be done 
at the level of amplitudes and generators deal usually with probabilities. 
As far as we know there is only one implementation based on the specific 
assumptions concerning amplitudes, and it applies only for a single string 
fragmentation processes f^. The most widely used procedure modelling the 
BE effect in the popular PYTHIA generator is based on the prescription 
for shifting the final state momenta to produce an enhancement at small 
momentum differences in the distributions of pairs of identical hadrons 
In this procedure one fits the parameters of the "input BE function" 

F{Q) = 1 + \exp{-Q^R^) (1) 

(assumed to have the same form as the standard parametrization of BE effect) 
to reproduce the experimentally observed effect. 

There is no simple relation between the values of input parameters A and 
R and the analogous parameters describing the experimental distribution. 
There is also no theoretical justification for this procedure and it should be 
regarded as a convenient parametrization, rather than the physical descrip- 
tion of the BE effect. 

The alternative approach is based on the formalism of Wigner functions 
0. One approximates the corrected distribution as a product of distribution 
without the BE effect and the weight function for which a definite prescription 



is given |[T^. This allows us to produce the distributions with the BE effect by 
generating the events without this effect and attaching to them the weights. 
To calculate these weights one must adopt several simplifying assumptions 
||10| (and hope they do not destroy the validity of the formulae). Finally one 
must assume the form of "two particle weight factor" and fit its parameters 
to describe correctly the data. 

Superficially, there is a marked similarity between these two approaches. 
In both cases the form of an "input function" is assumed and its parameters 
should be fitted to describe the data. However, there are also clear differences. 
Whereas the "momentum shifting method" has no theoretical justification, 
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it has two free parameters (plus a few hidden parameters defining the choice 
of momenta to be shifted, which affect the results rather mildly) and is quite 
easy to apply. Since all generated events are used, neither the multiplicity 
distributions nor, e.g., the decay channel probabilities in Z° decay are affected 
by shifting. There seems to be no need to change the values of the generator 
parameters fitted to the data before taking the BE effect into account. 

On the other hand, the weight method is quite well justified (granting 
that simplifying assumptions are not too rough), but there are many tech- 
nical problems with its use. Some of them have been solved: prohibitive 
increase of computational time with multiplicity may be avoided by a proper 
clustering procedure for final state momenta and the distortion of the 
multiplicity distribution may be removed by simple rescaling of weights de- 



pending on the event multiplicity [T^. Obviously, the weights may in prin 



ciple affect other distributions which were fitted to data without taking the 
BE effect into account. Thus the proper procedure would be to refit all the 
generator parameters comparing the weighted results with data. However, if 
the rescaling guarantees that average weight is equal one for each well defined 
class of events (e.g. in each Z° decay channel), the changes in distributions 
should be minor. 

Another notorious technical problem for weight methods is the instability 
of results due to the long tail of very high weight values. Usually it requires 
some arbitrary cut, but for sufficiently high number of generated events the 
effects of this cut are not very significant. Finally, there is a problem of 
selecting the particles, whose momenta are used to calculate weights and a 
problem of proper choice of "two particle weight factor" and its parametriza- 
tion (refiecting somehow the shape and the size of the production source). 

In this paper we discuss the solutions to the last two problems presenting 
the MC results for the BE effect in the hadronic decay of Z° and comparing 
them with some data. We consider only the distributions in the invariant 
four momentum difference 

Q' = -{pi-P2r. (2) 

The effect of anisotropy in various components of Q ||l3l was discussed else- 
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The following chapter contains the discussion of possible particle selec- 
tions and the infiuence of various MC parameters. The effects of different 
forms of two particle weight functions (considered already in the earlier pa- 
per H^) are presented in the third chapter. The last chapter presents the 
comparison with some data and conclusions. It should be stressed that we 
discuss only the standard prescription for weights justified by the Wigner 
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function formalism j^, |T^]. Other proposals ||T2|, |T^, should be regarded 
as viable versions of the weight method only if it is shown that they reproduce 
approximately the results obtained for this prescription. 



2 Particle selection and MC parameters 

Before discussing the details of the MC procedure implementing the BE effect 
we should decide which distribution will be used to present this effect. The 
standard quantity called "the BE ratio" is defined as 

^^^^^^ = 

where was defined above (0) and the numerator and denominator repre- 
sent the identical two-particle distribution with and without the BE effect, 
respectively. Obviously, this definition requires a more precise formulation 
of how we shall define the denominator. 

In the experimental definition of the BE ratio one uses often the distri- 
bution of unlike sign pion pairs but this requires cutting off the resonance 
effects. Thus recently it is preferred to use the pairs of identical pions from 
different events 

i^..«)^C-W).^^, (4) 
where the denominator is a convolution of single distributions 

Pi ® PiiQ) = I Pi{Pi)piiP2)5{Q^ + {Pi - P2?). (5) 



This choice of the denominator has other fiaws (i.e. it removes all correla- 
tions, and not only the BE effect). Therefore one uses often double ratios, 
dividing the experimental ratio by an analogous ratio of distributions from 
MC generator (without the BE effect) 

For the MC generated events the simplest choice is just to run MC without 
the procedure implementing the BE effect 

AIC,BE,^^ 

T3MC(r,\ — P"^ I7\ 
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Obviously if in the experimental investigation a double ratio is used, it seems 
more proper to calculate for comparison an analogous double ratio from MC 
events 



Fortunately the difference between R'^^{Q) and Rj^^{Q) is often in- 
significant. This is illustrated in Fig.l, where we show both ratios calculated 
for pion pairs from decays using the weight method with the Gaussian 
two-particle weight factor 



^-'2(^1,^2) = exp{ - ) (9) 



with a = 0.05 GeV^. 
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Fig.l. Comparison of the ratio of Q distributions (j^ with and without weights 
(squares) with double ratio of spectrum-to-background ratios with and without 
weights (stars). 

Here, as in all the later figures: 

• one million of events was generated by the PYTHIA 6.2 generator [0, 
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• the background distributions were constructed using pairs from differ- 
ent events; four million pairs of events were used for this purpose, 



• the BE ratios were normalized to approach smoothly the value of one 
at Q exceeding 1 GeV . 

For completeness, let us remind here that the two-particle weight factor 



is related to the full weight calculated for each event by a formula ||T0 



vr(pi,...,p„) =^n^2bi,pp(i)) (10) 

1=1 

where the sum extends over all permutations of n elements. More precisely, 
the event weight is a product of such sums calculated for all kinds of identical 
bosons registered by the detectors (in practice it is enough to count only 
positive and negative pions). 

This formula suggests that the a parameter in formula (|) is the only free 
parameter in the method. This would be, however, an oversimplification. 
The calculation of the full sum over permutations is practically impossible 
for the number of identical pions exceeding twenty [|l^. Thus we define 
the clusters of pions "close to each other" in the phase space and sum over 
permutations within clusters only. To make the results independent on the 
cluster definition we have to choose the value of the relevant parameter e much 
bigger than a. The details of this procedure have been described elsewhere 

Moreover, if the two particle weight factor is interpreted as a Fourier 
transform of the space-time distribution of pairs of pions, it seems justified 
to use a common shape of this factor only for pairs of "direct" pions. 

The decay products of long living particles and resonances are born far 
away from the original source and the corresponding two-particle weight fac- 
tor for pairs including these decay products would be close to the Dirac delta 
function, contributing negligibly to the final event weight. The same reason- 
ing was presented by Sjostrand who choose 20 MeV as a limit defining 
"long living" resonances and performed the momentum shift only for pions 
produced directly and the decay products of broader resonances. 

This suggests that we should calculate the event weight including in the 
sum only "direct" pions defined in a similar way. To avoid the changes of the 
original Monte Carlo procedures (which was the case for Sjostrand PYBOEI 
procedure, called internally from the generator before the decay of "long 
living" resonances and particles), we form a table of momenta for "direct" 
pions defined in various ways and use this table for the weight calculation. 
We found that the modifications of the original width limit of 20 MeV are 
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irrelevant as long as we do not include the u) decay products in the weight 
calculations. Including uj decay products enhances strongly the BE ratio, as 
shown in Fig. 2 for the decay. 

To justify the choice of the width limit let us require that the "direct" 
pion and the pion from ui decay have approximately the same momenta. The 
maximal momentum of a pion in the decay of uj at rest is about 300 MeV/ c, 
and the most likely value is of the order of 100 MeV/c. This allows to 
estimate that the distance between "birth points" of such pious is of the 
order of 10 fm and the corresponding width in momentum space should be 
about 20 MeV, smaller than the typical resolution. This suggests that the 
decay products of u (as well as the decay products of narrower resonances and 
other unstable particles) should not be taken into account when calculating 
weights. 
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Fig. 2. Comparison of the ratio of Q distributions with and without weights 
(0) for weights calculated excluding (squares) and including (diamonds) oj decay 
products. 

However, this argument has some flaws. First, the BE ratio is defined 
as a function of and not of the three-momentum squared (thus it reflects 
the space-time and not just the space structure of source). Second, the 
distribution of the decay length is not Gaussian. Thus we should not expect 
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a Gaussian shape of the weight factor. Finally, excluding the decay products 
of narrow resonances is a very rough procedure; a better solution would be 
to use different two-particle weight factors for different pairs of pions (direct- 
direct, direct-resonance and resonance-resonance). Let us add that all this 
estimate is classical and does not take into account the possible quantum 
effects. Thus we should not treat the choice of "direct" pions excluding the 
uj decay products as definite. In fact, the uncertainty of this choice seems 
to be the biggest uncertainty of the weight method. If necessary, it may be 
used to describe the BE effect if the observed values of the BE ratio at small 
Q'^ are high. 

The other free parameters of the PYTHIA generator may also infiuence 
the weights and the resulting BE ratio. An example of this effect is shown 
in Fig. 3 where we compare the results for default values of PYTHIA param- 
eters and for the values fitted to the L3 data [0. Let us stress that the 
choice of "direct" pions (excluding the uj decay products) and the value of 
a parameter are the same in both cases, but the results are visibly different. 
This is probably mainly due to the suppression of rj and rj' mesons for the 
L3 parameters. 
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Fig. 3. Comparison of the ratio of Q distributions with and without weights 
((J) for weights calculated using default PYTHIA parameters (squares) and L3 
parameters (circles). 
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Finally, let us check the dependence of the results on the value of the a pa- 
rameter. Until now we were using the value of 0.05 GeV^ which corresponds 
to the average Gaussian source size of the order of 1 fm. In Fig. 4 we compare 
the results (with L3 parameters) for a = 0.05 GeV'^ and a = 0.07 GeV"^. We 
see that by increasing a (which corresponds to a decreasing source size) we 
increase the width of "BE peak" and slightly increase its height. 
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Fig. 4. Comparison of the ratio of Q distributions with and without weights 
^ for L3 parameters with weights calculated using a = 0.05 GeV^ (circles) ana 
a = 0.07 Gey 2 (triangles). 

A notorious problem of the weight method is a possible distortion of var- 
ious distributions (fitted previously to data) by the introduction of weights. 
First such a distortion was observed for the multiplicity distribution where 
the probabilities of high multiplicities were enhanced by weights. This was 
cured by rescaling the weights with a factor CA" ||12[ where n is a charge 
particle multiplicity (measured in experiment). The values of parameters C 
and A are fitted to restore the original values of n and the original normal- 
ization. With this method the weights cause only a moderate increase of the 
dispersion of the multiplicity distribution. 

The weights infiuence also the single particle momentum distribution, 
reducing slightly the width, but these effects are not very significant. More 
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important is the change in the distribution of unhke sign pairs of pions, as 
shown in Fig. 5. A similar effect was observed for Sjostrand's implementation 
method of the BE effect. It should be noted, however, that by including the 
uj decay products in the weight calculation we increase the R ratio for unlike 
sign pion pairs by a few percent only, whereas the ratio for like sign pairs 
increased by about 50%, as shown in Fig. 2. Thus it is possible to describe 
a big BE effect without distorting seriously the distribution for unlike sign 
pairs. 




Fig. 5. Comparison of the ratio of Q distributions ^ with and without weights 
for unhke sign pairs of pions. Weights are calculated using L3 parameters with 
a = 0.07 GeV^ excluding uj decay products (triangles) and using default parameters 
with a = 0.05 GeV^ including u decay products (diamonds). 



3 Choice of the two-particle weight factor 

In the former section we used always the Gaussian two-particle weight factor 
d^). Obviously, there is no reason why all the space-time and momentum 
distributions should be described by such simple functions. However, if we 
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restrict ourselves to the monotonically decreasing weight factors normalized 
to one at = 0, it is easy to show that the Gaussian choice results in 
a curve which is the smoothest one and resembles the data best. This is 
demonstrated in Fig. 6 where we compare the results obtained for the default 
PYTHIA parameters for the Gaussian weight factor (|^) with a = 0.05 GeV^ 
and for two other choices of the weight factor: 



step-like 

with the same value of a 



ill) 



exponential 

with a = 0.03 GeV^. 
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Fig. 6. Comparison of the ratio of Q distributions ^ for default PYTHIA 
parameters and three different choices of two-particle weight factor: gaussian 
(crosses), step-like (solid line) and exponential (diamonds). 
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We see clearly that the shape of the obtained BE ratio reflects the shape 
of the two-particle weight factor. This may be written as an approximate 
relation 

Rbe{~{pi ~ 1 + 0-^2(^1,^2) (13) 

where the value of c depends on the shape of the weight factor and the 
selection of particles used in the weight calculations. 

Let us note that the Gaussian parametrization is unlikely to describe the 
data where the distribution of the space-time distance between the "birth 
points" of two pions is more complicated. This is the case for the four jet 
decay of W^W~ final states if two pions originate from two W-s. There it is 
unjustified to expect monotonically decreasing and normalized two-particle 
weight factors. However, for the Z° decay the parametrization (|^) seems to 
be the appropriate one. 

4 Comparison with data and conclusions 

We will not attempt here a detailed fit to any published data. There are 
many reasons for this reservation. First, as we have already mentioned, 
different experiments use different definitions of the reference sample in the 
denominator of the BE ratio. Thus the fit quality may depend on many 
factors not related to the procedure implementing the BE effect (e.g., the 
resonance effects and other correlations). Second, the published data include 
usually the acceptance corrections which are difficult to reconstruct in our 
calculations. In fact, the Monte Carlo parameters should be also fitted to the 
particular set of data before implementing the BE effect. As shown in Section 
2 there is a difference between the results obtained using default PYTHIA 
parameters and the parameters used by the L3 collaboration. 

Therefore we want to make only a semi-quantitative comparison between 
the results from our procedure and some high statistics data. To this purpose 
we use the recent L3 data shown as the reference sample in the paper devoted 
to the analysis of the WW decay ||T^. We compared them with various MC 
results shown in previous sections, rescaled with arbitrary constants to agree 
with data at > IGeV"^. We found that the modelled BE effect is too 
small compared with the data unless we include the u decay products for 
the weight calculation. We show the comparison in Fig. 7 for two choices of 
the a parameter in the weight factor (^. Normalization of both curves was 
adjusted to fit the data. We see that the data are qualitatively described by 
the PYTHIA MC with our implementation of the BE effect. One should not 
expect a good quantitative fit to the data for any single value of a; as already 
noted, one could use at least different values of this parameter (and, even 
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better, different sfiapes of W2) for pairs of pions of different origin. Then, 
however, the number of free parameters would increase making the success 
of the fitting procedure rather trivial. 
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Fig. 7. Comparison of the ratio of Q distributions for L3 parameters (Q) with 
a = 0.04 (soUd line) and a = 0.03 (broken Une) with the L3 data (stars). 

To summarize, we have discussed the freedom of the weight method for 
implementing the BE effect into Monte Carlo generators. We have shown that 
this freedom seems to be sufficient to describe the data. For pions coming 
from a single source which may be parametrized with a Gaussian distribution, 
there are three steps for choosing the weight method parameters: 

1. One should decide which ratio is used to display the BE effect and to 
calculate the same ratio from the MC with weights. It is preferred 
to use double ratios (cf. (|) and (||)) where both for data and MC 
one divides the chosen BE ratio by the same ratio calculated from MC 
without weights. 

2. One should choose the selection criterion for pions used to calculate 
weights. The typical choice corresponds to using direct pions and decay 
products of broad resonances, F > 20 MeV (as in Sjostrand's method). 
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3. One should select a proper value of the parameter a in (H). 
Technically, our algorithm contains four Fortran procedures: 

• LWBOEI, where for each event the "direct" pions are selected and their 
momenta are stored in the tables, 

• KLASKF, where pions of one sign are assigned to clusters, 

• PERCJE, where a weight factor from each cluster is calculated, 

• CLUSWAGI, where the full event weight is calculated as a product of 
weight factors from all clusters and all pion signs. 

All these procedures are available at request from us, together with a sample 
program calling the PYTHIA 6.2 generator and comparing the weighted and 
unweighted distributions for hadronic Z° decays. A modification of this pro- 
gram for other processes or other MC generators would be straightforward. 

One should also remember that after the introduction of weights one 
should rescale them by a CA" factor to restore the original normalization 
and average multiplicity. This, however, does not infiuence significantly the 
shapes of the BE ratios. 

To describe the process in which pions originate from two or more in- 
dependent sources (as the e"'"e~ — ^ W^W~ process with double hadronic 
decay of W-s) one needs a more elaborate procedure. Different forms of the 
W2 factor should be used for pairs coming from the same and from different 
sources. This will be discussed in detail elsewhere. 
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